Exact relaxation dynamics of a localized many-body state in the ID Bose gas 
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Through an exact method we numerically solve the time evolution of the density profile for an 
initially localized state in the one-dimensional bosons with repulsive short-range interactions. We 
show that a localized state with a density notch is constructed by superposing one-hole excitations. 
The initial density profile overlaps the plot of the squared amplitude of a dark soliton in the weak 
coupling regime. We observe the localized state collapsing into a flat profile in equilibrium for a 
large number of particles such as A*' = 1000. The relaxation time increases as the coupling constant 
decreases, which suggests the existence of off-diagonal long-range order. We show a recurrence 
phenomenon for a small number of particles such as. N = 20. 
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Localized wave packets play a fundamental role in 
quantum mechanics, particularly, in the study of dynam- 
ical behavior. For a massive particle in a free space it 
is possible to construct a localized state such as Gaus- 
sian wave packets by superposing single-particle eigen- 
functions, i.e., plane waves. When it evolves in time ac- 
cording to the Schrodinger equation, the wave packet col- 
lapses due to the uncertainty relations [l[ . For interacting 
many-body systems, on the other hand, it is in general 
a formidable task to pursue the unitary time evolution 
for a given initial state over a sufficiently long time with 
satisfactory accuracy 0] ■ It is even elusive how to define 
and construct an initially localized state of a physical 
meaning in a system with translation symmetry. 

The quantum dynamics of interacting many particles 
has recently attracted much interest associated with the 
question of "relaxation" and "thermalization" of isolated 
systems Q • In order to understand the concept of relax- 
ation von Neumann's ergodic theorem should be useful 
0-01 ■ Furthermore, relaxation dynamics after a sudden 
quench has been studied for integrablc systems Q. In 
the studies it is also essential to well define an initial 
state. Recent experiments in cold atoms Q enable us to 
study the isolated quantum dynamics of initially local- 
ized states because of the high isolation from the environ- 
ment and controllability of system parameters. Thus it 
is of importance to theoretically explore how to describe 
a localized state and to track its long-time behavior, in 
both integrable and nonintegrable systems. 

In this Letter, we show that a localized many-body 
state with a density notch is constructed by superpos- 
ing the Bethe ansatz eigenstates of a certain type for the 
ID interacting Bose gas with repulsive delta-function po- 
tentials. They are given by a series of one-hole excita- 
tions, which are called Lieb's type II excitations [l^l- In 
the weak coupling case, the initial density profile is quite 
similar to the graph of the squared amplitude of a dark 
soliton, which is a solution of the nonlinear Schrodinger 
equation for the classical scalar field. We numerically 



solve the exact time evolution of the expectation value of 
the density operator. Observing the movies of the density 
profile, we find that the system with N — 1000 shows re- 
laxation dynamics as if it were thcrmodynamically large, 
while that with = 20 shows finite-size effects such as 
recurrence phenomena. The relaxation time of the bro- 
ken symmetry state resembling a dark soliton increases in 
the weak coupling case. Thus, the localized wave packets 
are more stable as the constant c becomes smaller, which 
suggests the existence of off-diagonal long-range order. 

Let us consider the Hamiltonian of the ID interact- 
ing bosons with repulsive delta-function potentials, called 
the Lieb-Liniger (LL) model (lo| : 
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Here we assume the periodic boundary conditions of the 
system size L on the wavefunctions. Hereafter we con- 
sider the repulsive interaction: c > 0. The LL model 
is characterized by a single parameter 7 := c/n, where 
n — N/L is the density of the particles. We fix the par- 
ticle density as n = 1 and vary the coupling constant c. 
We employ a system of units with 2m = h = 1. The unit 
of time in our simulation is proportional to L~^. 

In terms of the canonical Bose field ilj{x,t), the 
LL model corresponds to "Hnls = /q dx[dxip^ dxi^ + 
cll^^^l''^^^l'^[^]. The second quantized Hamiltonian Hnls 
leads to the nonlinear Schrodinger equation for the 
canonical Bose field: idfiilj = —d'^'ip + ^dp"^ 'ip'ip . 

In the weak coupling limit (c — )■ 0), the one-particle re- 
duced density matrix is well approximated by the macro- 
scopic wavefunctions: pi{x,y) ~ '^*{x)'^{y), where 
^'(a;) is the eigenfunction corresponding to the maximum 
eigenvalue of pi{x,y). In fact, we can numerically show 
that the largest eigenvalue of pi is much larger than the 
other eigenvalues for small c, which suggests the exis- 
tence of off-diagonal long-range order. The density op- 
erator p{x) is defined by the diagonal elements of the 
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one-particle reduced density matrix: p{x) = pi{x,x). 

There is a well-known conjecture claiming that, in the 
weak coupling limit, the ID Bose gas, which is a quan- 
tum integrable system, should become a classical inte- 
grable system, which is described by the Gross-Pitaevskii 
equation, i.e., the nonlinear Schrodinger equation for the 
classical scalar field. It was addressed that the mode 
of classical dark solitons is identified with the type II 
excitations of the LL model through the coincidence of 
their dispersion relations in the weak coupling limit [lH . 
However, it is still nontrivial to show how the classical 
solitons are derived from the ID interacting bosons with 
c > even in the weak coupling limit. For the attractive 
case (c < 0), this was studied analytically [12[. Through 
numerical simulation, the quantum dynamics of dark soli- 
tons in the optical lattice system has been investigated 
for the Bose-Hubbard model 13| . The type II excitations 
are interpreted as quantum solitons in association with 
yrast states 

In the LL model, the Bcthc ansatz offers an exact 
eigenstate with an exact energy eigenvalue for a given 
set of quasimomcnta fci , fc2 , . . . , fcjv satisfying the Bethe 
ansatz equations for j ~ 1,2^ . . . , N: 
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Here Jj's are integers for odd and half-odd inte- 
gers for even N . We call them the Bethe quantum 
numbers. The total momentum P and energy eigen- 
value E are written in terms of the quasimomcnta as 
P = Y.U = X i:f=i I,. E = Ef.i kl If we spec- 
ify a set of the Bethe quantum numbers Ii < ■ ■ ■ < In, 
Bethe ansatz equations ([2|) have a unique real solution 
ki < ■ ■ ■ < kN [H]. In particular, the sequence of the 
Bethe quantum numbers of the ground state is given by 
Ij =■ -(TV l)/2 -h j for integers j with 1 < j < TV. 
The Bethe quantum numbers for low lying excitations 
are systematically derived by putting holes or particles 
in the ground-state sequence. 

Let us now construct an initial state with a localized 
density profile. In the type II branch, we denote by \P) 
the normalized Bethe eigenstate of total momentum P = 
2iTp/L for each integer p in the set {0, 1, . . . , — 1}. 
The Bethe quantum numbers of \P) are given by Ij = 
— {N + l)/2 + j for integers j with 1 < j < N — p and 
Ij = -{N + l)/2 + j + 1 for j with N - p+1 < j < N. 
For each integer q satisfying < q < N — 1, we define 
the coordinate state \X) of X = qL/N by the discrete 
Fourier transformation: 



\X) = 



1 
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p=0 



It turns out, remarkably, that the density profile of \X) 
has a density notch at the position X + L/2 as shown 
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FIG. 1: (Color online) Snapshots of the exact time evolution 
for the density profile with coupling constant c = 0.01 for 
iV = 1000, L = 1000 are shown (red solid lines). Also plotted 
in the top panel is the squared amphtude of a dark soliton 
(blue dashed line). 



in the top panels of Figs. [T] [21 and [H In the simula- 
tion we put g = 0, and the initial density profile p(x) is 
localized at x = L/2. The "temperature" of the state 
|X(i)) is estimated by equating the mean energy of the 
state with the internal energy derived from the thermo- 
dynamic Bethe ansatz [lij . although it is not in thermal 
equilibrium. 

The construction is analogous to the analysis of the 
attractive case (c < 0) in which the Fourier transforma- 
tion of the n-particle bound state |n, P) with momentum 
P gives a localized state \n,X) with the center of mass 
located at X, and the matrix element of the field opera- 
tor (n,X'|i/)(a;,t)|n-|-l,X) exactly corresponds to a bright 
soliton solution of the nonlinear Schrodinger equation for 
the classical scalar field with c < [l2| . 

Let us introduce \X(t)) = exp(— We evaluate 
the expectation value at time t of the density operator 
p[x) = [x^-^l^ix) for the initial state \X) as 



(X(t)|p(x)|X(t)) = l 



-2TTi{p-p')qlN 



p,p'—0 

X g»(p-p')x-»(iSp-iJp,)*(p'|p(o)|p), (4) 

where P = 2ttp/L and P' = 27rp'/L denote the total 
momenta of the normalized Bethe eigenstates \P) and 
|P'), respectively. 
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f = Q |N=L = 1000 
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FIG. 2: (Color online) Snapshots of the exact time evolution 
for the density profile with coupling constant c — 100 for 
N = 1000, L = 1000 are shown (red solid lines). Also plotted 
in the top panel is the squared amplitude of a dark soliton 
(blue dashed line). 



Because of quantum intcgrability, we numerically ob- 
tain all the energy eigenvalues Ep of |P)'s in the type 
II branch and follow the time evolution for quite a long 
time. 

We evaluate the form factor {P'\p{0)\P) in Eq. (g]) 
through Slavnov's formula [l7| together with the Gaudin- 
Korepin formula [3l for the norm of the Bethe eigenstate 
as 
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where the quasimomenta {fci, • • • , kN} and {k[, • • • , fc^} 
give the eigcnstates \P) and \P'), respectively. We use 
the abbreviations fcj.f :~ kj — kg and kjf := k'j — k'^. 
The kernel K{k) is defined by K{k) = 2cf{k^ + c^). The 
matrix G(fc) is the Gaudin matrix, whose (j, t) th element 
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1, 2, • • • , iV. The matrix elements of the by 
matrix U{k, k') are given by 



U [k, k')j,i = 2Sjilm 
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We thus numerically solve the exact time evolution of 
the density profile: Once we evaluate the form factors 
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FIG. 3: (Color online) Snapshots of the exact time evolution 
for the density profile with c = 100 for A*' = 20, L = 20 are 
shown. The solitary dip appears again at t = 33.0695. 



at t = we obtain the density profile at any late time 
t by taking only the sum of exponentials [l9[. Through 
Eq. ([5]) , the evaluation of the form factors of the Bethe 
eigenstates with N particles is reduced to that of the 
determinants of by iV matrices. 

The density profile of a wave packet for A^ ~ 1000 
initially has a localized shape without translation sym- 
metry, and then it relaxes into a flat profile in equilibrium 
through time development, as shown in Figs. [T] and [2] 
The relaxation time increases as the constant c decreases. 

In the weak coupling case of c = 0.01, the density pro- 
file is smooth and its time evolution is very slow as shown 
in Fig. [T] The density profile maintains a featureless 
curve through the time development and has a long co- 
herence in space suggesting the existence of off-diagonal 
long-range order. 

We observe that the initial density profile for c = 0.01 
overlaps the plot of the squared amplitude of a dark soli- 
ton solution of the nonlinear Schrodinger equation for the 
classical scalar field. This observation is fundamental for 
making an explicit connection between the type II exci- 
tations and the dark solitons. We shall discuss this in 
detail in a separate publication. 

In the strong coupling case of c = 100, the time evolu- 
tion is much faster. The initial density profile also shows 
small oscillating behavior at the shoulders of the density 
notch, which arc similar to the Friedel oscillations of the 
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FIG. 4: (Color online) The width of a density notch for the 
initial density profile and the healing length = \fcn with 
n = 1 are plotted against the coupling constant c. In the 
evaluation of the width of the density notch, we extrapolate 
the finite A'' data by 1/A'' expansion and estimate the values in 
the thermodynamic limit N ^ ca. The ratio between width 
and the healing length approaches a constant value in the 
weak coupling limit c — 0. 
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FIG. 5: (Color onhne) The Loschmidt echo F{t) = 
|(X(t)|X(0)}|^ is plotted against time t. The number of par- 
ticles are = 1000 (left panel) and iV = 20 (right panel). 
The particle density is N/L = 1 and the coupling constant is 
c — 100 for both panels. 



Fermi gas. The initial wave packet dynamically splits 
into many fragments, which results from the strong re- 
pulsive interaction. The initial density profile is rather 
different from the graph of the squared amplitude of a 
dark soliton. 

The width of the initial density notch is approximately 
proportional to the healing length ic for small c, as is the 
case for the classical dark soliton (see Fig. U]). For larger 
c, however, it seems that the width does not decrease in 
proportion to £c- This suggests the existence of fermionic 
hard cores. 

Let us calculate the Loschmidt echo, i.e., the fidelity. 
We define it by the overlap between the initial state 
|X(0)) and the time-evolved state \X{t)) at time t as 



F{t) ■.= \{x{t)\xm\'^^ 



N-1 
p=0 



(7) 



We plot the Loschmidt echo in Fig. [5] with fixed par- 
ticle density n = N/L = 1 in the strong coupling case of 
c = 100. For a large system size {N = 1000), as shown 
in the left panel of Fig. [51 it decays almost monotoni- 
cally with the short-time fluctuations being rather small. 



This behavior is common for a wide range of values of c 
with different scales of time. In the strong coupling limit 
c — > oo, by sending A^^ and i to oo with fixed density n = 
N/L we derive analytical forms of the Loschmidt echo: 



Fit) 



S{nV^)^]/i2TrnH), where C(a 



and S{x) are given by the Fresnel integrals defined by 
C(x) := /(j^ cos(7rs^/2)ds and S{x) := sin(7rs^/2)ds, 
respectively. As a short-time expansion we have F{t) = 
1 — Air^n'^t'^ /45 + 0{t^), and as a long-time expansion 
we have F{t) = {AwnH)-^ + 0{l/t^). The plot of nu- 
merical results for N = 1000 with c = 100 (left panel of 
Fig. [5]) almost overlaps that of the analytical expression 
of the Fresnel integrals. It suggests that the system of 
N = 1000 already shows thermodynamic behavior. 

For a small system size {N = 20) the Loschmidt echo 
is shown in the right panel of Fig. [5j Here, we observe 
large fluctuations due to the finite-size effect. 

Let us define a long-time average E by E[/(t)] := 
ImiT^oo dtf{t)/T. Using E[e*"*] = for cj 7^ and 
E[e*'^*] = 1 for cj = 0, we obtain E[F(t)] - 1/7V for any 
value of c. This is also confirmed by our simulation. 

A sharp peak at i ~ 33 in the Loschmidt echo for 
= 20 shows the signal of a recurrence phenomenon. 
In fact, the localized wave packet is revived at this time 
after a "tentative" relaxation as shown in Fig. |31 

In conclusion we have shown that the superposition of 
the Bethe eigenstates of the type II excitations leads to a 
quantum many-body state with a localized density pro- 
file. It perfectly overlaps the plot of the squared ampli- 
tude of a dark soliton in the weak coupling regime. By 
means of Slavnov's formula together with the Gaudin- 
Korcpin formula, we have numerically solved the exact 
time evolution of the density profile in both the strong 
and weak coupling cases for a large number of particles 
such as A^ = 1000 over a very long period of time. For 
a sufficiently large system size (A^ = 1000), we observe 
that a density profile with broken translation symmetry 
relaxes into a flat profile in equilibrium through the time 
development. For a small system size (A^ = 20), the local- 
ized wave packet is revived after a tentative relaxation. 
These observations suggest that the method presented 
in this Letter should be fundamental for exploring ex- 
act approaches to the quantum dynamics of many-body 
systems. 
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